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ABSTRACT 

Since z ~ 1, the stellar mass density locked in low mass groups and clusters has grown by a factor 
of ~ 8. Here we make the first statistical measurements of the stellar mass content of low mass 
X-ray groups at 0.5 < z < 1, enabling the calibration of stellar-to-halo mass scales for wide-held 
optical and infrared surveys. Groups are selected from combined Chandra and XMM — Newton 
X-ray observations in the Chandra Deep Field South (CDFS). These ultra-deep observations allow 
us to identify bona fide low mass groups at high redshift and enable measurements of their total 
halo masses. We compute aggregate stellar masses for these halos using galaxies from the Carnegie- 
Spitzer-IMACS (CSI) spectroscopic redshift survey. Stars comprise ~ 3 — 4% of the total mass of 
group halos with masses 10 128 < M 2 oo/Mq < 10 13 ' 5 (about the mass of Fornax and l/50th the mass 
of Virgo). Complementing our sample with higher mass halos at these redshifts, we find that the 
stellar-to-halo mass ratio decreases toward higher halo masses, consistent with other work in the local 
and high redshift universe. The observed scatter about the stellar-halo mass relation is a ~ 0.25 dex, 
which is relatively small and suggests that total group stellar mass can serve as a rough proxy for halo 
mass. We find no evidence for any significant evolution in the stellar-halo mass relation since z < 1. 
Quantifying the stellar content in groups since this epoch is critical given that hierarchical assembly 
leads to such halos growing in number density and hosting increasing shares of quiescent galaxies. 


1. INTRODUCTION 

The buildup of stellar mass in dark matter halos has 
been measured in the local universe with a variety of 
methods, indirectly w i th abundance match ing techniques 
(jBehroozi et all 120131 iMoster et all 1201 3li and more di¬ 
rectly with meas urements of stellar and total masses of 
individual halos dKravtsov et all 12014 ). These stellar 
mass - halo mass relations provid e stro ng constra ints on 
models of galaxy formation fe.g.. iGenel et alJl2014l) . 

While much of the attention has been focussed on the 
stellar content of L* halos, recent work has highlighted 
the relevance of the higher mass scales represented by 
groups and clusters, whe re various feedback p rocesses 
are predicted to operate. iKravtsov et al.l (120141) for in¬ 
stance show discrepancies between stellar-to-halo mass 
ratios computed for clusters from abundance matching 
and from direct observations. Meanwhile, at z < 1, 

I Williams et al.l (120121) measure a dramatic increase in the 
number density of galaxy groups over an epoch when star 
formation activity is g reatly diminished, bot h in relative 
intensity (i.e., SSFR, iFumagalli et alJ 12012ft and in the 
number of galaxies that display active levels of star for¬ 
mation at a fixed mass. Given that these massive halos at 
high redshift represent sites where galax ies are q uenched 
at levels above that in the field (IPatel et al.ll2009t 120111) . 
an analysis of environmental trends provides insight into 
the mechanisms that contribute toward the growth of the 
red sequence. Directly measuring the mass scales of these 
halos in which galaxies undergo quenching will provide 
stronger constraints on the relevant mechanisms than in- 
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direct methods for inferring halo membership and halo 
mass scales. 

The two relevant measurements for directly probing 
group mass scales are the halo mass and the total group 
stellar mass. At high redshift, reliable detection of 
galaxy groups as well as halo mass estimates require 
deep X-ray imaging. We therefore utilize the ultra deep 
Chandra+XMM New ton g roup catalog in the CDFS 
of iFinoguenov et all (|2014l ). To find galaxies within 
these groups, we require a survey with sufficient red¬ 
shift precision, which is uniform, and which has well- 
characterized completeness corrections to relatively low 
stell ar masses. Th e Ca rnegie-Spitzer-IMACS (CSI) sur¬ 
vey dKelson et al.ll2014l) meets these criteria and is used 
in this work. 

We assume a cosmology with H 0 = 70 km s _1 Mpc -1 , 
Dm = 0.30, and = 0.70. Stellar masses are based on 
a IKroupal (f2f)0lh IMF. 

2. DATA AND ANALYSIS 

2.1. Chandra and XMM-Newton X-ray Imaging 

We utilize combined Chandr a and XMM-Newton X- 
ray observations in the CDFS (|Finoguenov et aD 120141) 
to select groups at high redshift. These observations rep¬ 
resent the deepest X-ray data to date. For the analy¬ 
sis p resented here, we cull th e list of 46 X-ray groups 
from IFinoguenov et al.l (12014 ) to those with redshifts in 
the range 0.5 < z < 1 and flags of 1 or 2, indicat¬ 
ing well defined centroids and redshift determinations 
from followup spectroscopy. In addition, we include only 
those groups with CSI coverage in the ECDFS, out of 
4.8 sq degs in the SWIRE CDFS field. Only 23 groups 
remain in the ca t alog. Total masses, M 2 oo, given in 
IFinoguenov et all (|2014l ) are estimated from the Lx — M 
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Figure 1 . Example X-ray detected groups in the CDFS. The color images are constructed from HST/WFC3 J 125 and H 160 bands. The 
centers of the groups are marked with a green plus sign and the virial radii, R 200 , are indicated by the large dashed green circles. Galaxies 
that are characterized as members of the group are circled in blue, while nearby field galaxies are circled in red and orange based on their 
redshift offset from the group. Clear overdensities exist in the vicinity of the identified group centers. 


relation in lLeauthaud et all (l2010f) . which was calibrated 
with weak lensing based mass estimates in the COSMOS 
field. Though the lowest halo masses were derived by ex - 
trapolating the latter relation, iFinoguenov et al. ; 201 1: 
confirmed the reliability of the halo mass estimates from 
stacked weak lensing and clustering analyses. i ?200 val¬ 
ues follow from these mass estimates. 

2.2. CSI Spectroscopy and Completeness 

We utilize a redshift catalog in the CDFS constructed 
from new CSI observations. The handlin g of the data 
is sim ilar to that in the XMM-LSS field (jKelson et all 
MB), though our deep DECarn photometry provides for 
more uniform completeness as a function of magnitude. 
Briefly, CSI targets were selected in the Spitzer IRAC 
3.6 /jm band with magnitudes brighter than 21 AB mag 
and observed with the Uniform Dispersion Prism (UDP) 
on IMACS (|Dressler et alJ[2011fl at Magellan. Redshifts 
were measured by fitting the available spectrophotome¬ 
try (spectrum + ugrizJHK) with SEDs generated from 
stellar population synthesis mo dels . In comparisons with 
high-resolution redshifts from I Cooper et all (120121) we 
find a biweiglrt scatter of <r z /( 1 + z) ~ 0.008 for both 
red and blue galaxies. At 2 = 1, the stellar mass com¬ 
pleteness limit is M ~ 10 10 3 M 0 . 

We note that group numbers 79 and 68 do not contain 
CSI galaxies above the mass cut of M > 10 10 3 M 0 (some 
galaxies may have been missed due to slit collisions and 
other sources of incompleteness), though both contain 
galaxies below this limit. These are both the lowest mass 
groups in the sample with halo mass estimates of Mh ~ 
10 12 ’ 8 M 0 . Also, number 63 had a large uncertainty on 
the stellar mass estimate (see below). We exclude these 
groups from our analysis leaving a total sample of 20 
X-ray groups in the CDFS. 

The spectroscopic completeness for CSI in the CDFS 


is ~ 50 — 60% over 17 < i < 25.5 AB mag and a broad 
range of color. To accurately account for galaxies that 
were not observed with the UDP or which did not have 
well constrained redshifts, galaxies are assigned weights 
according to their inverse completeness in bins of ob¬ 
served magnitude ([3.6]), obs erved c olor (i — [3.6]), and 
local source density (see, e.g.. Kelson et al.ll2014T) . 

Even though the CDFS has been extensively studied by 
other surveys, with numerous high -resolu tion red shifts 
available in the literature (e.g.. ICooner et ahll2012f ). we 
refrain from combining the multiple available datasets. 
The different selection functions in other work would 
make it difficult to produce a homogeneous dataset from 
the aggregate of all the available surveys, despite poten¬ 
tially boosting the level of completeness. The depth and 
homogeneity of the CSI data simplify the analysis, and 
the relatively uniform ~ 50 — 60% completeness is high 
enough to robustly characterize the stellar mass content 
of these X-ray groups. Lastly, while other works have 
used photometric catalogs to assemble or cali b rate g alaxy 
samples in groups in th e CDFS ('e.g.- IGiodini et afll2009t 
lErfanianfar et ahll2014f) . our sample is exclusively spec¬ 
troscopic. 

2.3. Group Membership and Aggregate Stellar Mass 
Calculation and Uncertainty 

For each X-ray group, we assign CSI galaxies above the 
mass limit within the projected virial radii, R 2001 and 
within |Az|/(l + z) < 0.02. This large redshift window 
encompasses both members and contaminants (see below 
for correction). We sum the stellar masses of the galaxies 
down to M s tars(> 10 10 - 3 M 0 ), multiplying by the com¬ 
pleteness weights. On spatial scales comparable to those 
of the groups at 0.5 < z < 1 (e.g., D ~ 1'), the com¬ 
pleteness correction accurately reproduces galaxy counts 
to within ~ 20%, an uncertainty that is factored into the 
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Figure 2. Selection properties of the X-ray group sample, (a) Rest-frame X-ray luminosity (0.1-2.4 keV) versus redshift for groups from 
our primary sample in the deep CDFS field (black) and the comp ariso n sample in the shallower COSMOS field (gray). Filled circles indicate 
flag=l while open symbols indicate flag=2 (see text in Section 1 3.1 1) . (h) M 200 versus redshift for the same sample. The CDFS sample 
reaches lower halo masses at 2 > 0.5, probing below M < 10 13 ' 3 Mq. 


total error budget on the aggregate stellar mass measure¬ 
ments. 

To account for galaxies below our stellar mass com¬ 
pleteness limit, we correct these group stellar masses as¬ 
suming a Schechter function with a = — 1 and log M* = 
10.9. Integrating down to M = 10 9 M 0 we find a cor¬ 
rection factor of 1.27. Changing a and M* within rea¬ 
sonable ranges yields an uncertainty of ~ 30% in this 
correction, which is added to the error budget. 

We derive a contamination correction for every group 
by subtracting off the stellar mass contribution from fore¬ 
ground and background galaxies. We measured the stel¬ 
lar mass within circular apertures of radius R 200 at spa¬ 
tial locations that did not overlap with other groups near 
%-oup- We compute the aggregate stellar mass of galax¬ 
ies at |Az|/(1 + Zgroup) < 0.02 in these apertures and 
take the mean and standard deviation to represent the 
background and its uncertainty (typically ~ 18%). 

After measuring the stellar mass within I? 200 > we cor- 
rect the value to that expected within r 2 oo assuming a 
spherically symmetric NF W profi le (Nav arro et al.lfl996f) 
seen in projection (e.g., IGiodini et all 120091 . see their 
Equation 7). We assume an NFW profile with concentra¬ 
tion of c = 4 resulting in a multiplicative correction fac¬ 
tor of ~ 0.78. The range of typical concentr ation values 
assoc iated with high mass halos (c ~ 3 — 5, iDuffv et~ahl 
I2008D impacts the correction factor by only ~ 3% and 
therefore does not significantly affect our results or con¬ 
clusions. 

Lastly, there is an additional source of uncertainty 
due to the statistical uncertainties in the individual 
mass measurements. The typical statistical uncertainty, 
based on repeat UDP observations and spectrophotomet- 
ric SED fitting, is ~ 0.1 dex per galaxy. This does not 
include contributions due to systematic effects such as 
different stellar population synthesis models and IMF 
variations. 

The total uncertainty in the aggregate stellar mass 
measurement for a typical group in our sample is ~ 


0.2 dex. 

Figure [T] shows two example low-mass groups in the 
CDFS. Color images were constructe d usin g HST CAN- 
DELS WFC3 J 125 and Lfifio imaging (IGrogin et al.ll20lit 
IKoekemoer et a l. 2011). Within the virial radii (dashed 
green circles), Figured] shows clear overdensities of group 
members (blue circles) relative to field galaxies within 
the same redshift slice (red and orange circles), indicat¬ 
ing that we are able to reliably place galaxies into groups 
with our data and compute robust aggregate stellar mass 
totals. 

2.4. COSMOS Comparison Sample 

We compare our results for low-mass groups at 0.5 < 
z < 1 to low-mass groups at 0.1 < z < 0.5 and to high- 
mass groups at 0.5 < z < 1, both derived from the wider 
field, but shallow er X-ray data from C OSMOS (for a de¬ 
tailed study, see IGiodini et all I2009D . We standardize 
the group stellar masses in COSMOS with ours in the 
CDFS by measuring them in the same manner as de¬ 
scribed in the previous section and using the photomet¬ 
ric redshifts and stellar masses from lGeorge et al.1(l2011l ). 
While this approach differs from summing redshift likeli¬ 
hood functions within the reds hift windows of the groups, 
as carried out in IGeorge et all (]2011l) . there is very good 
agreement between the two methods with a scatter in the 
aggregate stellar masses of only ~ 0.06 dex and a neg¬ 
ligible offset. We convert the COSMOS stellar masses 
from Chabrier to Kroupa IMF by adding 0.05 dex. In 
total, the COSMOS sample consists of 125 groups, 69 at 
0.5 < 2 < 1 and 56 at 0.1 < z < 0.5 with halo masses 
spanning 10 12 8 < M 2 oo/Mq < 10 14 3 and with X-ray 
flag=l or flag=2 and no other flags. 

3. RESULTS 

3.1. Galaxy Group Sample Properties 

Figure [2jz. shows the X-ray luminosity of X-ray groups 
as a function of redshift. Redshifts were determined 
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Figure 3. The stellar mass-halo mass relation (M s tars versus M 200 ) for high mass halos at 0.5 < z < 1. Stellar mass measurements 
represent the aggregate mass in stars of all galaxies in the halo, (a ) Individual X-ray detecte d groups in the CDFS and COSMOS are in 
black and gray, respectively. Clusters from GCLASS are in green (Ivan der Burg et aI.II20l4h . Measurements and uncertainties from the 
three datasets have been standardized when possible. Dotted lines across the figure indicate ratios of stellar-to-halo mass of 0.01 and 0.1 for 
reference. A biweight fit to the CDFS+COSMOS data (blue line) indicates a slope of 0.84 ± 0.10 and an observed scatter of a = 0.25 dex. 
Fitting all of the data, including GCLASS (green line), the slope becomes 0.72 ±0.07 and the observed scatter cr = 0.25 dex. The less-than- 
unity slope indicates that fewer stars are formed/assembled in higher mass halos. The relatively small observed scatter implies that total 
stellar mass serves as a good proxy for the total halo mass in this regime. (6) Median stellar to halo mass ratio for different bins in halo 
mass. Vertical error bars represent bootstrapped uncertainties. Black horizontal error bars indicate the bin range while the thicker gray 
error bars represent the 25-75th percentile values for the halo masses in each subsample. At low group masses of 10 12 ' 8 < M/Mq < 10 13-5 
probed mostly by the data in the CDFS, stars make up 3 — 4% of the total mass of the halos. The orange curve indicates the z 0 relation 
from abundance matching. 


from followup spectroscopy of red seq ue nce c andi¬ 
date galaxies as d escrib ed in iGeorge et al.1 (|201lH and 
iFinoguenov et all (12014 submitted). The filled symbols 
represent X-ray groups with centroids determined from 
the X-ray image (flag=l), while the open symbols repre¬ 
sent centroids determined from the positions of the red 
sequence galaxies (flag=2). The median redshift for the 
CDFS sample is z me a ~ 0.7. Figure© shows the COS¬ 
MOS and CDFS group masses, M 2 oo, as a function of 
redshift. The error bars represent the statistical uncer¬ 
tainty and do not include the contribution due to the 
scatter in the L\ — M relation (~ 0.2 dex). As a re¬ 
sult of the lower flux limits, CDFS reaches significantly 
lower group masses at high redshift compared to COS¬ 
MOS. In particular, halos at z > 0.5 with mass below 
M 2 oo < 10 13 - 3 Mq are almost exclusively detected in the 
CDFS. 

3.2. The Stellar Mass - Halo Mass Relation at 
0.5 < z <\ 

The aggregate stellar mass within X-ray groups at 
0.5 < 2 < 1 is compared to their total mass in Fig¬ 
ure [U In Figure ©, the CDFS groups are shown us¬ 
ing the black points. The COSMOS groups are shown 
in light gray. Owing to its depth, the CDFS sample 
provides the only leverage on M sta rs/Af 2 oo for low mass 
groups (M 2 oo < 2 x 10 13 Mq). 

The blue line is a biweight fit to the CDFS and COS¬ 
MOS data of the form 

log M stars = /? + a log M 200 (1) 

where a = 0.84±0.10 and /? = 0.47±1.33 (see Table[l]). 


The slope is consisten t with the value of a = 0.81 ± 
0.11 found in Giodini et all (2009) for their higher mass 
COSMOS group sample spanning 0.1 < 2 < 1. The 
implied stellar fraction for the COSMOS groups (~ 2% 
at M?no ~ IQ 13,5 M g) ) is a lso consistent with the results 
of iLeauthaud et al.1 ( 2012 :. The green data points in 
Figure [3] represent measur ements for clust ers from the 
GCLASS survey (|van der Burg et all I2014H . The green 
line is a fit to all three datasets and results in a = 0.72 ± 
0.07 and /3 = 2.20 ± 0.97. The inclusion of the cluster 
data results in a shallow slope which is consistent with 
the general trend found at high halo masses. 

The observed biweight scatter about the relation in 
Equation Q] is a = 0.25 dex, and is comparable for both 
the CDFS and COSMOS samples. The similar scatter is 
reassuring and implies that our completeness corrections 
are robust and unbiased. Furthermore, this observational 
scatter is remarkably small. When factoring in the un¬ 
certainties on the aggregate stellar mass measurements, 
the intrinsic scatter must be quite low. Determination 
of the intrinsic scatter requires accurate estimates of the 
observational errors. Because they total ~ 0.2 dex, we 
are unable to make a precise determination of the intrin¬ 
sic scatter. In any case, an upper limit for the intrinsic 
scatter in the group stellar mass versus halo mass rela¬ 
tion of ~ 0.25 dex reflects a relatively tight correlation 
suggesting that the former serves as a rough proxy for 
the l atter. This has also been noted by other authors 
(e.g., Ivan der Burg et alll2014f) . 

Figure© shows the median stellar to halo mass for dif¬ 
ferent bins in halo mass. The bulk of the low mass CDFS 
group sample lies in the two lowest mass bins spanning 
10 12 ’ 8 < M 2 oo/M 0 < 10 13 5 . In these halos, which are 
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Figure 4. Ratio of group stellar mass to halo mass versus redshift 
for two bins of halo mass. Low mass halos (10 12 8 < M/M@ < 
10 13,3 ) are shown in blue and high mass in red (10 13,3 < M/Mq < 
10 14 ). Results from abundance matching at 2 : ~ 0 J Kravtsov et al.l 
[2014) are shown as open squares. We do not find any significant 
evolution in the ratio of stellar to halo mass for groups at 2 < 1. 


Table 1 

Best fit parameters to stellar-halo mass relation 


Sample Fitted 

a 

P 

&obs 

(dex) 

N 

CDFS, COSMOS 

0.84 ±0.10 

0.47 ± 1.33 

0.25 

89 

CDFS, COSMOS, GCLASS 

0.72 ±0.07 

2.20 ±0.97 

0.25 

93 


Note. — The fit coefficients a and /3 are indicated in Equation[T] 
log Mstars = /3 a log M 200 • The observed biweight scatter about 
this relation is cr 0 b s and the number of halos included in the fit, TV. 
Fits apply only to the mass range probed by the data presented in 
this work. 

10-20 times more massive than the Milky Way’s halo, the 
stellar mass constitutes 3 — 4% of the total mass budget. 

3.3. No Significant Evolution in the Group Stellar to 
Halo Mass Ratio since z ~ 1 

Here we examine the redshift evolution, if any, of the 
stellar mass - halo mass relation for groups. The or¬ 
ange curve in Figure^ indicates the stellar-to-halo mass 
r atio from abundance matching, empirically derived by 
iKravtsov et al.l 1 2014 Figure 12) at 2 ~ 0. Given the er¬ 
ror bars and any additional uncertainty due to systematic 
differences between the two measurement methods, there 
does not appear to be any significant evolution at z < 1 
in the stellar-to-halo mass relation for M 2 oo > 10 13 M 0 . 

We can further examine and confirm the lack of evo¬ 
lution by using additional measurements at lower red- 
shifts (0 < z < 0.5) that were made in the same man¬ 
ner as our own. Figure 2] shows the redshift evolu¬ 
tion for the median stellar-to-halo mass ratio for two 
halo mass bins. Uncertainties were computed by boot¬ 
strapping. Note that in the lowest mass bin, 10 12 ' 8 < 
M 200 /M 0 < 10 13 ' 3 (blue), the entire sample at z > 0.5 is 
comprised of groups in the CDFS where the X-ray imag¬ 
ing is of sufficient depth to detect these systems and stel¬ 
lar masses have been measured with CSI. The low red¬ 
shift (z < 0.5) sample is compri sed exclusively of COS¬ 
MOS X-ray groups (see Section Gl!. The open squares 


indicate the results a t z ~ 0 from abundance matching 
(|Kravtsov et al.lf2014f) and the associated error bars the 
range of median values spanned by halos with masses in 
the 25-75th percentile of the two sub-samples. Again, we 
do not find any significant evolution in the stellar-to-halo 
mass ratio in the group regime at z < 1. 

4. DISCUSSION 

We have measured the stellar mass - halo mass relation 
for low mass groups at 0.5 < z < 1, complementin g 
other work at these reds hifts (e.g.. iGiodini et aLll2009t 
Ivan der Burg et all 12014 ). For the low mass groups in 
our CDFS sample, we find that about ~ 3 — 4% of the 
total mass is contained in stars, less than that of L* 
haloes but more than that of clusters, as is also the trend 
in the local universe. The lack of (or slow) evolution in 
the stellar fraction of such halos is relevant given that the 
number density of these grou ps is observed to increase 
dramatically at z < 1 (William s et, al.ir2012li indicating 
the growing significance of such environments. 

Late-time formation of massive halos is expected from 
hierarchical growth. Massive halos in the mass range 
M 200 > 10 13 M© are predicted to comprise ~ 32% of 
the total mass of all halos above M 200 > 10 11 ' 6 M© (ap¬ 
proximately the halo mass scale that would host central 
galaxi es at our stellar mass limit) at z ~ 1 (ba sed on HM- 
Fcalc. iMurrav et afll2013t [Tinker et al.ll2008l b By z ~ 0, 
this percentage rises substantially to ~ 52%. For groups 
in the mass range 10 13 < M 200 /M Q < 10 14 , the percent¬ 
age increases from ~ 29% at z ~ 1 to ~ 36% at 2 ~ 0. 
Clusters at M 2 00 > 10 14 M 0 clearly contribute toward 
the mass growth at z < 1 bu t even t hese struc tures are 
built, in part, from groups (lBerrier et al.1 1200 9). 

The increase in number densities of galaxy groups, 
and the associated increase in the stellar mass density 
within such halos, has implications for the buildup of 
quiescent galaxies (QGs) since z ~ 1. From z ~ 1 
to z ~ 0 the quiescent fraction (QF) increases from 
~ 44% to ~ 59% for galaxies with M > 10 10 3 M 0 
dKelson et al.l 120141 : iMuzzin et al.l 120131) . Down to this 
mass limit, QFs are h igher in rich environments com¬ 
pared to the field (]Dresslerill980f ) with QFs of ~ 80% and 
~ 35% in groups/clusters and the field, respectively (e.g., 
iPatel et ad1l201lh . As a consequence, the rise of groups 
over the past ~ 8 Gyr contributes substantially to the 
evolving QF over that time. By assessing the star for¬ 
mation activity and quiescence of galaxies in groups over 
the full 15 degs 2 survey field of CSI, we will accurately 
constrain the impact of large scale structure growth on 
the buildup of QGs and the decline in global SFRD since 
z ~ 1. 

We thank Alexis Finoguenov for making his CDFS X- 
ray catalog available. 
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